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^ Abstract 

C. 

(D | Within the framework of FRW cosmology, we study the QCD modified ghost 

scalar field models of dark energy in the presence of both interaction and viscos- 
ity. For a spatially non-flat FRW universe containing modified ghost dark energy 
(MGDE) and dark matter, we obtain the equation of state of MGDE, the decel- 
eration parameter as well as a differential equation governing the MGDE density 



Of 



parameter. Furthermore, we reconstruct both the dynamics and potentials of the 
quintessence, tachyon, K-essence and dilaton scalar field DE models according to 
the evolution of the MGDE density. 
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1 Veneziano MGDE paradigm 



The Veneziano MGDE density is given by [1] 

PD = aH + f3H 2 , (1) 

where a and (3 are two constants, roughly of order of Aq CD and Aq CD , where Aqcd ~ 100 
MeV is the QCD mass scale and H is the Hubble parameter. The original version of the 
Veneziano GDE, i.e. po = has been already studied in the literature [2, 3, 4, 5, 6, 7, 
8, 9, 10]. Note that in Eq. (1), the term proportional to H 2 becomes important during 
the early evolution of the universe and behaves like the early DE [1]. 

In a spatially non-flat FRW universe filled with MGDE density p^ and DM energy 
density p m , the first Friedmann equation within the framework of Einstein gravity reads 

where the scalar curvature k = 0,1,-1 denote a flat, closed and open FRW universe, 
respectively. Also M P = (8-RGy 1 ^ 2 is the reduced Planck mass. 
Using the fractional energy densities 



n _ Pm _ Pm n _PD__ PD r, _ « ( o\ 

m ~ Pcr ~ 3MJ,# 2 ' p cr _ 3Mj,if 2 ' k ~ am* ' [6) 



the Friedmann equation (2) can be rewritten as 

i + n k = n D + n m . (4) 

Substituting Eq. (1) into the middle relation of Eq. (3) gives 



a 



where 



7 = 1 - (6) 



3M 



Inserting this into the last relation of Eq. (3) yields 

n x J ^^-)\ (7) 

("D + 7 - 1) V a ) 

where we take a = 1 for the present time and the subscript "0" denotes the present value 
of a quantity. 

Following the observational evidences we extend our study to the case in which the 
MGDE has a bulk viscosity property [11] and also interact with DM [12]. In the presence 
of bulk viscosity and interaction, the energy densities of MGDE and DM do not conserve 
separately and continuity equations take the forms 

p D + 3H{l + u D )p D = 9H 2 £-Q, (8) 



2 



p m + 3Hp m = Q, 



(9) 



where up = Pd/pd is the equation of state (EoS) parameter of MGDE. Also £ = eH~ 1 p£ ) 
is the bulk viscosity coefficient with the viscosity constant e and Q = 3b 2 H(p D + p m ) is 
the interaction term with the coupling constant b 2 . 

Taking time derivative of Eq. (2) and using Eqs. (3), (4), (8) and (9) one can get 



H 

IP 



1 + ^ + (u D - 3e)n D 



(10) 



Taking time derivative of Eq. (1), using (10), and substituting the obtained result into 
Eq. (8) gives the EoS parameter of the interacting viscous MGDE as 



1 



u D = 



n D -i-i 

For the deceleration parameter 



i _ 9i + 1 ~ 1 

3 Q D 



1 + - ) + 2b 2 1 1 + ^ 



n_ k 

3 



D 



+ 3e. 



H 



(12) 



replacing Eq. (11) into (10) gives 

7 + 2(Q D - 1) - Q k + 3b 2 (l + Q k ) 

Taking the derivative of Eq. (36) with respect to In a and using Eqs. (10) and (11) gives 



(13) 



D 



din a 



n D - 1 -i t 



D -i-Y + b2 ( 1 + ^ 



(14) 



It is interesting to note that Eqs. (13) and (14) show that the deceleration and MGDE 
density parameters are independent of viscosity constant. 

We solve the differential equation governing the MGDE density parameter, Eq. (14), 
numerically. In Fig. 1, variation of the MGDE density parameter versus redshift 
z— ^ — 1 for different coupling constants b 2 is plotted. Figure 1 shows that: i) For a given 
b 2 , Qd increases during history of the universe, ii) At early and late times, Qd increases 
and decreases, respectively, with increasing b 2 . 

In Fig. 2 we plot the evolutionary behavior of the deceleration parameter, Eq. (13), for 
different b 2 . Figure 2 presents that: i) the universe transitions from a matter dominated 
epoch at early times to the de Sitter phase, i.e. q = —1, in the future, as expected. The 
result for b 2 = is in agreement with that obtained by [1]. ii) For b 2 = 0.0, 0.02 and 
0.04 at z — 0.75, 0.91 and 1.10, respectively, we have a cosmic deceleration q > to 
acceleration q < transition which is compatible with the observations [13]. iii) For a 
given z, although q decreases with increasing b 2 , in the far future q becomes independent 
of6 2 . 

The variation of the EoS parameter of interacting viscous MGDE, Eq. (11), for dif- 
ferent coupling b 2 and viscosity e constants is plotted in Figs. 3, 4 and 5. Figures 3 and 4 
illustrate that: i) in the absence of viscosity (e = 0), ujd(z) varies from the quintessence 
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phases (lod > — 1) to the phantom regime (u>d < — 1)- ii) In the absence of interac- 
tion {b 2 = 0), wo(z) varies from uo > —1 to oj d = —1, which is similar to the freezing 
quintessence model [14]. iii) For a given z, up decreases and increases with increasing b 2 
and e, respectively. Note that the result for b 2 = e = is in accordance with that obtained 
by [1]. Figures 5 shows that the interaction and viscosity have opposite effects on uup. 

In the next sections, we establish a correspondence between the MGDE and various 
scalar field models by identifying their respective energy densities and equations of state 
and then reconstruct both the dynamics and potential of the field. 



2 Quintessence MGDE 



The quintessence scalar field energy density and pressure are given by [15] 

1 



PQ 



Pq 



2 

1 ; 
2 



2 - v{4>). 



The quintessence EoS parameter takes the form 

_ Pq _ 4> 2 ~ 2V(<P) 
UQ Pq <P 2 + 2V(4>y 



(15) 
(16) 

(17) 



Identifying Eq. (17) with the EoS parameter of interacting viscous MGDE (11), ooq = wo, 
and also equating Eq. (15) with (1), pq = p D , one can get 



V{<t>) 



[l-u D )p D , 



4> 2 = (l + u D )p D . 



(19) 



Inserting Eqs. (1) and (11) into the above equations, one can get the quintessence poten- 
tial and kinetic energy as 



v{4>) 



a 2 Q 



D 



+ f -(i±)(l 



3 



n D - 7 - i 



3e 



,(20) 



2 



a 2 n 



D 



3M|(Od + 7 -1) ; 



+ (i + f ) + 2&2 (^) 



n D - 7 - 1 

Integrating this gives the modified quintessence scalar field as 

"^-7-f + (l + f)+2^(ig 



+ 3e 



(21) 



0(a) - 0(a;) = V3M P [ a 



n B - 7 - 1 



+ 3e 



da 



(22) 
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The above integral cannot be taken analytically. But with the help of Eq. (7) and numer- 
ical solution of the differential equation (14) one can obtain the evolutionary behavior of 
the quintessence MGDE scalar field. In Figs. 6, 7 and 8 we plot the variation of scalar 
field, Eq. (22), versus redshift. Figures 6 and 7 present that: i) for a given b 2 or e, 
increases during history of the universe, ii) For a given z, decreases and increases 
with increasing b 2 and e, respectively. Figure 6 shows that for b 2 = 0.0, 0.02 and 0.04 at 
z < —0.97, —0.57 and —0.41, respectively, becomes pure imaginary, i.e. 2 < 0, and 
does not show itself in Fig. 6. For 2 < the MGDE scalar field behaves as a phantom- 
type scalar field [16]. Figure 8 clarifies that interaction and viscosity have opposite effects 
on 0. For b 2 = e, the effect of viscosity is more than the interaction. 

In Figs. 9, 10 and 11, the variation of quintessence MGDE potential, Eq. (20), versus 
redshift is presented for different b 2 and e. Figures 9 and 10 illustrate that: i) for a given b 2 
or e, the quintessence MGDE potential V(4>) decreases during history of the universe, ii) 
For a given z, V(4>) increases and decreases with increasing b 2 and e, respectively. Figure 
11 shows that in the presence of both interaction and viscosity, although for b 2 = e at 
early times the effect of viscosity is more than the interaction, at late times they neutralize 
the effect of each other. 



3 Tachyon MGDE 



The energy density and pressure of the tachyon scalar filed are as follows [17] 



V(<f>) 

\/l-0 2 ' 



Pt = 



PT = -V(0)Vl-0 2 . 

The EoS parameter of the tachyon filed reads 



ojt 



- 1. 



Pt 



From Eqs. (1) and (23), p D = p T gives the kinetic energy term 

^ 2 _ ^-7-f + jr(i + f)+26 2 (^) ^ 

Q D - 7 - 1 

Also using Eqs. (11) and (25), ujd = wt yields the tachyon potential 



a 2 Q 



D 



3M 2 {Q D + 1 - l) 2 



3 



2b 2 (2±0*) 



-3e 



1/2 



(23) 
(24) 
(25) 

(26) 
(27) 



Integrating Eq. (26) with respect to a and using (36) gives the evolutionary form of the 
modified ghost tachyon scalar field as 

0(a)-0(a,) = ^ / (n D + 7 -l) 

OL Jen 



X 



Vn-l-f + g(l + f)+26 2 



+ 3e 



-• (28) 
a 
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In Figs. 12, 13 and 14 we plot the variation of 0, Eq. (28), versus redshift for different b 2 
and e. Also the evolution of tachyon MGDE potential, Eq. (27), is shown in Figs. 15, 16 
and 17. Figures 12 to 17 clarify that the scalar field and potential of the tachyon MGDE 
model behave like the quintessence one (see Figs. 6 to 11). 



4 K-essence MGDE 

For the K-essence scalar field [18, 19], the pressure 

PK(<p,x) = fm-x+x 2 ), 

and energy density 



give the EoS parameter 



M0,X) = /(0)(~X + 3X 2 ), 



Pk X ~ 1 
Pk 3x - 1 



(29) 
(30) 
(31) 



where x — 2 /2 is the K-essence kinetic energy. Following Eqs. (11) and (31), ud = Wff 
gives 

2 + 7 -f + ^(l + f)+26 2 (i±^)-^ + 3e(fi D - 7 -l) 

X = 



4 + 7 - &k + g(3 + n k ) + 66 2 (igt) - Q D + 9e(fi D - 7 - 1) ' 
Integrating this with respect to a gives the scalar field of the K-essence MGDE as 

f(fi D + 7 -l) 

J Obi 

'7-1 



(32) 



0(a) - 0(oj) = 



3MJ, 



x 



4 + 2 7 



1 + f ) + 46 2 (i±^) - 2fi D + 6e(Q D - 7 - 1) 



4 + 7 - fi fc + (g) (3 + Q k ) + 66 2 (i±0t) - Q D + 9e(fi D - 7 



-■(33) 



The evolutionary behavior of this field for different b 2 and e is plotted in Figs. 18, 19 and 
20. Figures present that the K-essence MGDE scalar field increases during history of the 
universe. 



5 Dilaton MGDE 

For the dilaton scalar field [20], the pressure and energy density take the forms 

p D = -X + ce A V, (34) 

PD = -X + 3ce A V- (35) 

Here c and A are two positive constants and x — 2 /2 is the dilaton kinetic energy. 
Therefore, the dilaton EoS parameter reads 

Pd -1 + ce^x ( „ R , 

UD = y D = -l + Sce^x (36) 

6 



Identifying Eq. (11) with (36) reduces to 



ce^x 



2 + 7 ~f + g(l + f)+26 2 (l0)-^ + 3e(^- 7 -l) 
4 + 7 - tt k + 2^(3 + Q k ) + 66 2 (i0) - fi D + 9e(fi D - 7 - 1) " 



(37) 



Using x = 2 /2 one can take the integral of above equation with respect to a. The result 
yields 

2 



0(a) 



A 



x 



+ 5 / (fi D + 7 -l) 



4 + 2 7 



+ 2 



2: 



r)(i + 



46 2 



2/1+fli 



2fi D + 6e(fi D - 7 -l) 



4 + 7 - fi fc + (g) (3 + fi fc ) + 662 (i±Ot) _ fi D + 9e (Q D - 7 - 1) 

3M^A 



da 
a 



,(38) 



where 5 = 2^775- I* 1 Figs- 21, 22 and 23, variation of the dilaton MGDE scalar field, Eq. 
(38), is illustrated for different b 2 and e. Figures clear that the scalar field of the dilation 
MGDE model behaves like the K-essence one (see Figs. 18, 19 and 20). 



6 Conclusions 

Here we studied the Veneziano MGDE model in the framework of Einstein's gravity. We 
considered a spatially non-flat FRW universe filled with interacting viscous MGDE and 
DM. We derived a differential equation governing the evolution of the MGDE density 
parameter and solved it numerically. We also obtained the EoS parameter of the in- 
teracting viscous MGDE and the deceleration parameter of the universe. Furthermore, 
using a correspondence between the interacting viscous MGDE and quintessence, tachyon, 
K-essence and dilaton scalar field models of DE we reconstructed the dynamics <fi and po- 
tentials V{4>) of the aforementioned scalar field models according to the evolution of the 
MGDE density. Our numerical results show the following. 

(i) The MGDE density parameter for a given coupling constant b 2 , increases when 
the time increases. The evolutionary behavior of Qd is independent of viscosity. 

(ii) The variation of the deceleration parameter q shows that the universe transitions 
from an early matter dominant epoch to the de Sitter era in the future, as expected. Also 
q like Qn does not depend on viscosity. 

(iii) The EoS parameter un of the MGDE model in the absence of viscosity (e = 0), 
varies from the quintessence phase to the phantom regime. Whereas in the absence of 
interaction (6 2 = 0), it behaves like the freezing quintessence model. The interaction and 
viscosity have opposite effects oiiud. 

(iv) The quintessence and tachyon MGDE scalar fields for a given b 2 or e, increase 
during history of the universe. For a given redshift, they decrease and increase with 
increasing b 2 and e, respectively. The potentials of the quintessence and tachyon MGDE 
models for a given b 2 or e, decrease with increasing the time. The interaction and viscosity 
have opposite effects on both the scalar fields and potentials of the aforementioned models. 

(v) The K-essence and dilaton MGDE scalar fields for a given b 2 or e increase during 
history of the universe. 
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Figure 1: MGDE density parameter, Eq. (14), versus redshift for different coupling 
constants b 2 . Auxiliary parameters are f2 fco = 0.01, Q Do = 0.76 and 7 = 1.105 [1]. 
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Figure 2: Deceleration parameter, Eq. (13), versus redshift for different coupling constants 
b 2 . Auxiliary parameters as in Fig. 1 
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Figure 3: EoS parameter of MGDE, Eq. (11), versus redshift for different coupling 
constants b 2 with e = 0. Auxiliary parameters as in Fig. 1. 
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Figure 4: Same as Fig. 3 for different viscosity constants e with 6 2 = 0. 
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Figure 5: Same as Fig. 3 for different coupling b 2 and viscosity e constants. 
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Figure 6: Quintessence MGDE scalar field, Eq. (22), versus redshift for different coupling 
constants b 2 with e = 0. Auxiliary parameters are f2 fco = 0.01, £l Do = 0.76, 7 = 1.105 [1] 
and 0(1) = 0. 
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Figure 7: Same as Fig. 6 for different viscosity constants e with b 2 = 0. 
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Figure 8: Same as Fig. 6 for different coupling b 2 and viscosity e constants. 
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Figure 9: Quintessence MGDE potential, Eq. (20), versus redshift for different coupling 
constants b 2 with e = 0. Auxiliary parameter as in Fig. 6. 
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Figure 10: Same as Fig. 9 for different viscosity constants e with b 2 = 0. Auxiliary 
parameters as in Fig. 6. 
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Figure 11: Same as Fig. 9 for different coupling b 2 and viscosity e constants. Auxiliary 
parameters as in Fig. 6. 
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Tachyon MGDE (e=0) 
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Figure 12: Tachyon MGDE scalar field, Eq. (28), versus redshift for different coupling 
constants b 2 with e = 0. Auxiliary parameters as in Fig. 6. 
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Figure 13: Same as Fig. 12 for different viscosity constants e with b 2 = 0. 
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Figure 14: Same as Fig. 12 for different coupling b 2 and viscosity e constants. 
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Tachyon MGDE (6=0) 




Figure 15: Tachyon MGDE potential, Eq. (27), versus redshift for different coupling 
constants b 2 with e = 0. Auxiliary parameter as in Fig. 6. 
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Figure 16: Same as Fig. 15 for different viscosity constants e with b 2 = 0. 
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Figure 17: Same as Fig. 15 for different coupling b 2 and viscosity e constants. 
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K-essence MGDE (e=0) 
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Figure 18: K-essence MGDE scalar field, Eq. (33), versus redshift for different coupling 
constants b 2 with e = 0. Auxiliary parameters as in Fig. 6. 
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Figure 19: Same as Fig. 18 for different viscosity constants e with b 2 = 0. 
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Figure 20: Same as Fig. 18 for different coupling b 2 and viscosity e constants. 
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Dilaton MGDE (6=0) 
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Figure 21: Dilaton MGDE scalar field, Eq. (38), versus redshift for different coupling 
constants b 2 with e = 0. Auxiliary parameters are f2fc = 0.01, flo = 0.76, 7 = 1.105 [1], 

#1) = and* = = 1. 
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Figure 22: Same as Fig. 21 for different viscosity constants e with b 2 = 0. 
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Figure 23: Same as Fig. 21 for different coupling b 2 and viscosity e constants. 
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